function H = H1_matrix(h)
    Hx1_matrixH1 = h(1);  Hy1_matrixH1 = h(2);  Hz1_matrixH1 = h(3);
    H = [
        Hx1_matrixH1,   0, -Hz1_matrixH1, Hy1_matrixH1;
          0, Hx1_matrixH1,  Hy1_matrixH1, Hz1_matrixH1;
       -Hz1_matrixH1, Hy1_matrixH1, -Hx1_matrixH1,   0;
        Hy1_matrixH1, Hz1_matrixH1,    0,-Hx1_matrixH1];
end